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Abstract 

We extend the calculation of the Unruh effect to the universality classes of quantum vacua obeying topolog- 
ically protected invariance under anisotropic scaling r — > br , t — > b z t. Two situations are considered. The 
first one is related to the accelerated detector which detects the electron - hole pairs. The second one is 
related to the system in external electric field, when the electron - hole pairs are created due to the Schwinger 
process. As distinct from the Unruh effect in relativistic systems (where z = I) the calculated radiation 
is not thermal, but has properties of systems in the vicinity of quantum criticality. The vacuum obeying 
anisotropic scaling can be realized, in particular, in multilayer graphene with the rhombohedral stacking. 
Opportunities of the experimental realization of Unruh effect in this situation are discussed. 



1. Introduction 

Relativistic physics can be considered as emergent phenomenon, which occurs in the low-energy limit 
in some universality classes of quantum vacua with topologically protected nodes in the energy spectrum. 
The general topological approach to classification of gapless vacua, which is based on connection to the 
if-theory, has been developed by Hofava [l|. Hofava demonstrated that the reason for the emergence of 
Lorentz invariance is the topological theorem, called the Atiyah-Bott-Shapiro construction. 

The Lorentz invariance emerges only for the particular classes of topological vacua, where nodes in energy 
spectrum are characterized by the elementary topological charges, N — ±1. For the higher values of the 
winding number, the Lorentz invariance can be supported by symmetry between the fermionic species. This 
symmetry makes the vacuum to be effectively equivalent to system of fermions, each having elementary 
topological charges N — ±1 , see Q. If such symmetry is absent but the splitting of the nodes in the energy 
spectrum into the elementary nodes is still forbidden, then in the quantum vacua with higher values of 
topological number | N\ > 1 the low-energy fermionic spectrum obeys the invariance under anisotropic scale 
transformation. An example of the anisotropic scale transformation is r — > br , t — >• b z t with 1, i.e. one 
has different scale transformation for space and time and thus contrary to the relativistic systems with z = 1 
space and time cannot be united into space-time. Systems with anisotropic scaling in 3+1 and 2+1 vacua 
are discussed in Ref. Q and in Refs. 0, 0, [jj H, 0, [§], respectively. In a simple case of 2+1 systems, such 
as multilayered graphene 0, the scaling parameter z coincides with topological charge of the node in the 
energy spectrum, z = \N\, in other words, the anisotropic scaling is topologically protected. 

Physics obeying the topologically protected anisotropic scaling is fundamentally different from the rela- 
tivistic physics. It emerges in the vacua, which belong to different topological classes, and thus cannot overlap 
with the relativistic physics. The transition between different topological classes cannot occur adiabatically. 
That is why the transfer from the vacuum with emergent relativistic physics to the anisotropic vacuum with 
a given scaling parameter z ^ 1 occurs only by discontinuous jump - by quantum phase transition. The 
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same refers to the quantum transitions between anisotropic vacua with different z. Thus, the topologically 
protected quantum vacua with anisotropic scaling open the new natural area for application of quantum 
field theory. For example, while traditional relativistic QED was considered as the function of the space 
dimension D , the general QED must be considered as function of two parameters, D and N 8] . The 
QED with the anisotropic scaling can be realized in the multilayered graphene, and it is essentially richer 
than the relativistic QED with z = 1. The modified Heisenberg-Euler action describing the polarization 
of the anisotropic quantum vacuum with z ^ 1 has been discussed in Refs. @, EL 12|- In particular, this 
effective action has different scaling laws for electric and magnetic fields. 

The anisotropic scaling can be useful also to quantum gravity. For example, the vacua with topologically 
protected anisotropic scaling can host the Hofava gravity, which is based just on the idea of anisotropic 



scaling 13|, UM |l5(. The advantage of the Hofava extension of general relativity is that for z = 3 the 



quantum gravity can be ultraviolet complete. 

The natural question is how the other quantum phenomena, such as Hawking radiation 16[ and Unruh 
effect [17| . which are intrinsic to the quantum vacuum, are modified compared to the relativistic case. The 
related phenomenon - Schwinger pair creation in electric field in anisotropic vacuua - has been considered 



in Refs. |8J, [ill, [12J. In the relativistic vacuum, the Bckenstein conjecture of the generalized second law 



of the black hole thermodynamics [l8[ and Hawking conjecture of the black hole radiance [l6| indicated 
the close connection between general relativity (GR) and thermodynamics of the quantum vacuum. This 
was followed by the Unruh conjecture [l7[ that the detector moving in the Lorentz invariant vacuum with 
constant acceleration views the quantum vacuum as a thermal background. The latter indicated the close 
connection between the thermodynamics and the vacuum in special relativity. Then the thermal behavior 
of the quantum vacuum has been extended to the expanding Universe [l9j |. At the same time, it was found 
that the vacuum of Lorentz invariant quantum electrodynamics (QED) in an exactly constant electric field 
E can be also understood as a thermal background characterized by the temperature parameter [20j ] (for 
recent papers see [ill]). 

All these phenomena occurring in the relativistic quantum vacuum are interrelated and have similar 
properties. In particular, the power spectrum of the vacuum noise (or the detector-response function) seen 
by a uniformly accelerated observer in flat spacetimes depends on space dimension D, and for even D it 
exhibits the phenomenon of the apparent inversion of statistics (23j | . A similar phenomenon of the inversion 
of spin-statistics relation occurs for the vacuum polarization in a constant electric field [20L [25j j . As another 
example, it is argued that the same back reaction effect of particle creation influences the vacuum of QED 
in electric field and the quantum vacuum in the de Sitter space-time [26L This imp lies the instability of the 
de Sitter vacuum towards decay, which remains the issue of controversy [27 , 28 , 2t| 3(J 3ll : Extension of the 
stability analysis to de Sitter vacuum with different z 7^ 1 has been considered in Ref. [321 ] . 

Here we discuss the modification of the thermal properties of the quantum vacuum as compared to 
the Lorentz invariant vacua. We consider the particle production in electric field and the Unruh effect 
in the anisotropic scaling system, such as multilayered graphene. In general anisotropic case the energy 
distribution of electron-hole pairs in electric field and the distribution detected by the accelerated observer 
do not look as thermal. Nevertheless the distributions depend on the same dimensionless quantity E/T, where 
T is some effective parameter of dimension of energy, which is determined by electric field or acceleration 
cor resp on ding ly and by the integer parameter z. This indicates the new type of quantum criticality (33L EmI . 
I35L l36l I37I l38j emerging in the anisotropic vacua with z ^ 1. We will also discuss the opportunities to study 
these effects in the multilayer graphene experimentally. 



2. Unruh temperature in multilayer graphene 

In this section we extend the notion of Unruh temperature Tjj to the field theory of multilayer graphene 
(with ABC stacking). Namely, we introduce it as the typical energy of processes that take place in the 
presence of external force F acting on the quasiparticles. (In the case of external electric field £ we have F — 
e£.) We consider classical motion of quasiparticles, introduce the generalized acceleration that characterizes 
this motion, and relate it to the Unruh temperature. We give qualitative estimate of the "typical" energy in 
graphene in the presence of external force. This estimate is considered as the definition of Unruh temperature 
and will be given in Eq. (|7|. It will be shown in the following sections that this temperature enters the 
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considered distributions over energy u> of the electron - hole pairs in a number of situations. Namely these 
distributions depend on the ratio ui /Tjj. 



2.1. Classical motion of quasiparticles 

We deal with the two-component spinors placed in the external electric field directed along the x-axis. 
We introduce the external force F acting on the fermion excitation adding the potential energy Fx to the 
one-particle Hamiltonian. Then the Hamiltonian for low-energy fermionic quasiparticles propagating in the 
system with anisotropic scaling z = J has the form 

(0 v(p x - ipy) \ 

( \ J \ (1) 

V{p x +ipy) 0. I 

Here v — vp is the Fermi velocity for the case of monolayer (J = 1) and v — l/2m for the case of bilayer 
(J = 2), where m ~ 0.028m e [39[ and m e is the free-electron mass. In general case v oc t^ J v F , where tj_ is 
the interlayer hopping parameter. Here and further in the text we work in the system of units with h = 1. 
Schrodinger equation has the usual form 

idt$> = . (2) 
We consider the fermion excitations with fixed value of the transverse momentum p y : ty(t, x, y) — e tPyV ip(t, x). 
The conserved transverse momentum p y plays the role of mass of the one-dimensional fermion (multiplied 
by vp instead of the speed of light). We denote p y = M = Mvp. The corresponding energy spectrum 

is E — ±v(p x 2 + M 2 ) = ±tj^ v F (p x 2 + A4 2 v F ) . Since p y plays the role of mass, it is natural to 

identify the generalized acceleration as a — F/M — vpF/p y . To justify this choice, let us apply classical 
approximation to ip(t,x) with the corresponding classical Hamiltonian 

/ \ J/2 

H c i{p xi x) = Fx + v[p x 2 +p ) . (3) 



Hamilton equations of motion give classical trajectory that corresponds to i(0) = 0: 

x{t)=vF J - 1 (t 2 + V \) J/2 - vF' 1 - 1 ^ + x(0) , a=^-=v F -. (4) 
V a z / a J M. p y 

This is the generalization of the relativistic J = 1 case, where Eq.(j4|) describes the hyperbolic motion with 
correctly determined linear acceleration a. The particular case p y = corresponds to the trajectory 

x{t)^vF J - 1 \t\ J + x{Q) (5) 

2.2. Unruh temperature 

As it was mentioned above we define the Unruh temperature Tjj as the "typical" energy of processes 
taking place in the system under the action of the external force. We estimate this energy semiclassically 
using the above mentioned classical trajectory for the particular case p y = 0. Hamiltonian equations relate 
classical momentum to time as p c i — Ft. In Eq. ([3]) both terms of the Hamiltonian are typically of the same 
order. So, we can estimate the energy scale as Tjj ~ vp J cl ~ vt J F J . Here typical values of time are related 
to energy by the uncertainty relation tT\j ~ 1. Thus we get 

Tv-LfA^ 1 (6) 



For J — 1 we shall typically use the definition Tjj = (v F Fy /(2tt). For J > 2 we shall define T v = 

(yF J ^J J+1 . For the completeness let us rewrite the expression for the Unruh temperature at J > 2 in the 

usual system of units (with the dimensional Fermi velocity for the single-layer graphene vf and the Plank 
constant fi restored): 



Tu =t± 



with t± « 0.4 eV, and v F sa c/300, where c is the speed of light. Here it is used that v = t\_ J v J F . 
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2.3. Generalized acceleration 

The identification of Eq. ((4]) for the generalized acceleration does not work for the case of Eq. ([5]), i.e. for 
P y = 0. This is the case that was used above for the estimate of the "typical" energy Tjj. Unruh temperature 
in (quasi) relativistic models is related to the acceleration as Ty — a/(2 / KVp). This allows us to estimate 
the "typical" acceleration in the presence of the external force F for the multilayer graphene. For J > 1 we 
define the acceleration as equal to the Unruh temperature of Eq. ([7]) multiplied by the Fermi velocity vf- 

The same result is obtained when we define the generalized acceleration as F/M e f f , where the effective 
mass M e s is played by vacuum fluctuations of the p x projection of momentum, (Ip^D- The latter can be found 
from the Hamiltonian (|3]) at p y — 0. The contribution of two terms to the energy of vacuum fluctuations is 
of the same order, while the fluctuations (\x\) and (\p x \) are related by the Heisenberg uncertainty relation: 
(\x\) ~ l/(\ Px \). Then one has F/<|p x |> = v{\p x \) J which gives M eff ~ (\p x \)/v F = (F j Vp an d the 
generalized acceleration a = F/M e ff. 

In the special case J = 2 the equation a = F/M c g is applicable also for massive case M ^ 0. The reason 
for that is that the trajectory Q for J = 2 

Ft 2 

x(t) = vFt 2 + x(0) = — + x(0) , (8) 
2 m 

does not depend on M and thus the generalized acceleration should be the same as for M — 0. The tra- 
jectory ([8]) is similar to the trajectory of a Galilean particle with mass m, and at first glance the properly 
determined acceleration is the Galilean acceleration a — F/m. However, the considered particles are essen- 
tially non-Galilean, since they have positive and negative branches of energy spectrum E = ±v[p x 2 + p y 2 ^j , 
which correspond to the gapless energy spectrum with parabolic touching. That is why instead of Galilean 
acceleraton a = F/m, the characteristic acceleration is determined by the quantum fluctuation of p x on the 
trajectory, as in the case of M = 0, i.e. a = F/M e g ~ %f (1/mF) 1 ' 3 . 

So, in what follows the generalized acceleration in 2+1 anisotropic scaling systems can be defined as 

a = v F F/\p y \ for J^2 , Py^O, (9) 

a = vfF^ — J = v F T v for J = 2 , or for J > 2 , p y = . (10) 



3. Accelerated observer interacting with quantum fluctuations 

In this section we consider the effect of the vacuum fluctuations on the observer accelerated by the 
constant external force in the 2+1 quantum vacuum with anisotropic scaling. As an example, we shall 
use the multilayer graphene with rhombohcdral (ABC...) stacking [9|, in which the topologically protected 
scaling exponent z coincides with the number of layers, z = J 0, 

3.1. Trajectory of the detector, the interaction between the detector and the fermionic excitations, and the 
click rate 

In principle, one can use any trajectory of the detector. In this section we shall consider for J > 2 the 
trajectory given by Eq. ([5]). For the special case J — 1 we consider different trajectory because Eq. ([5]) at 
J = 1 would lead to the singularity. 

In this section we assume that the detector interacts with the fermion excitations via the term 

y interaction ~ ^ + O (11) 

with some operators Oa- Different choices of these operators correspond to different practical realizations of 
Unruh detector. In real graphene various experimental situations may be suggested corresponding to different 
practical realizations of the Unruh detector. In subsection !3.7l we shall briefly consider one of the realizations, 
- the Raman scattering from the spot moving along the graphene sheet. In this case rough approximation 
gives Vinteraction ~ ifi + cr a cr b ip (a, b — 1,2) if the scattered light is not polarized and its direction is orthogonal 
to the graphene sheet. It is worth mentioning that in this particular situation as well as in the other possible 
experiments the microscopic description of the Unruh detector is rather complicated. However, the properties 
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of the spectra become independent on the practical realization, and, in particular, on the form of Vinteraction 
for the energies much larger than Tjj, when the semiclassical approximation works. Therefore, expressions 
obtained semiclassically in subsection 13.61 are universal. Nevertheless, for the completeness we consider 
several cases, when the problem can be solved exactly. The considered cases correspond to idealized Unruh 
detectors corresponding to different trajectories and different interaction terms Eq. (|f f [) . 

The response function responsible for the transition of the detector to the excited state is (Eq. (3) of 



F t (u>) = 2Re / dfe-^W&t-f) (12) 



In Eq. (TT21 we encounter the Wightman function W(ti,t2) of two currents i/j + OaiP- However, for t\ > t 2 
it is equal to the T-ordered two-particle Green function. We represent it as 

W(h,t 2 ) = {Ti/j+ih, x[h])0 A ip{tx, x[h])ij + {t 2l x[t 2 })0+^{t 2 , x[t 2 ])) 

~ TyO A e(ti-ta,a:[ti]-a![t2])Oi0(t2-ii,rs[t2]-a;[ii]), h > t 2 (13) 

where Q is the single-fermion Green function. It is worth mentioning that the expression for the click rate 
written in this form resembles the conventional expression for the cross - section of the deep inelastic processes 
like e + + e~ — > hadrons. Both quantities are related to the imaginary parts of the corresponding polarization 
operators. In fact, the derivations of these expressions are similar. In Eq. (1131) we should consider carefully 
the limit t\ — t 2 — >• (see discussion at the end of this subsection). 

The summation over A is assumed in Eq. (fT"3"|) . Below we concern several particular choices of Oa ■ In 
pure relativistic 2 + 1 theory there are the following two natural choices: Vinteraction ~ ">pip = 4> +a ' 3 ' t l , an d 

Vinteraction ~ ^(T^lj) = 1p+ <7 3 'd^lp , (1=1,2,3. 

In non - relativistic models of graphene the natural choices are: Interaction ~ = tp + o~ 3 ipi ^interaction ~ 
4> + 4>, Vinteraction ~ i> + a "^ , a =1,2, and Vinteraction ~ ip + a a cr b , a,b = 1,2. The latter case corresponds to 
Raman scattering from the moving spot for the nonpolarized light incoming (and outgoing) within the beam 
orthogonal to the graphene plane (see subsection 13 . T|) . In the relativistic case in 4D the Unruh detectors 
interacting with the electron- hole pairs were considered in |47| . 

The response function is expressed as follows: 

F M~Re / dte- iut TrO A g(-t,-x[t])0+g{t,x[t}) (14) 
Jo 

We need to substitute into this expression the particular trajectory of the detector and the particular 
operators Oa corresponding to the interaction of electron - hole pairs with the internal degrees of freedom of 
the detector. The integral in Eq. ([M)) may be divergent at t — > 0. Let us consider the polarization operator 



n(t, a;) = iTr A Q{~t, -x)0\ Q(t, x) (15) 

p +oo 



The expression Q(u),x) — Im J Q °° dtexpy—iuitj TI(t,x) gives the probability that the detector at rest (x = 

const) absorbs (to > 0) or emits (u < 0) the electron - hole pair. Correspondingly, Q(ui,x) should vanish 
for lj > because there are no free electron - hole pairs in vacuum. We consider the analytical continuation 
of H(t, x) to negative values of t, and extend integration to the whole real axis of t. Now we already must 
remember, that we deal with the Wightman functions. According to the general properties of Wightman 
functions, U(—t,x) — — U*(t,x). We close the contour in the lower half of the complex plane. The resulting 
integral vanishes only if there are no poles within this contour. That is, n(i, x) is analytical in the lower half 
of the complex t - plane. 

If the singularities appear at t € R, we should go around them from the bottom. This gives the rules of 
going around poles. These rules are equivalent to the existence of infinitely small imaginary contribution to 
t. We are to substitute everywhere H(t — ie, x) instead of H(t, x). Thus we shift in Eq. (fT4|) the integration 
to t — it: 

/>oo 

F (uj)~Re dte- iut TrO A g(-t + ie,-x[t])0^g(t-ie,x[t]) (16) 



This is the standard way of the regularization of Wightman functions (see, for example, [22j, |24|). We adopt 



it throughout the text of the present paper when the nonrelativistic results are concerned. However, for the 



true relativistic case considered in section 13.41 there is a certain complication [45j, |4 
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3.2. The two point Green function 

The two point Green function can be expressed as 



Q(t,x)~ J ded 2 p e -^ t+ ^ x g(e,p) (17) 
Here 0(e,p) is the fermion Green function in momentum space: 



G ^P> = w~T7 = R[p] + rTTT^tP] 

e — H \p\ e — v\p\ J a A 



1 



(w) 




In order to go around the poles correctly in the integral over e we consider the Wick rotated Green functions. 
We introduce the Euclidean time z: t = —iz. Integration over the Euclidean frequency u (e = — iuj) gives 
for z > 0: 

duj -rr^- = 2m [ 6 n " I 19 

, w-iu|p| J CT 3 V / 

One can see, that for z > the Green function describes the propagation of the excitation with positive 
energy (electron). As a result 



G{iz,x) = -m J d 2 pe iM - vzM " 




Gn G12 
G\i Gn 



At z < we have 



Gn(iz,x) = — 2it 2 i J pdp J a (\p\x) exp(—vz\p\ J ) 

Gi 2 (iz,x) = -2t?i J pdpi J J.j(\p\x)cxp(-vz\p\ J ) (20) 

roo e iuz / Q Q \ 

^ 0,-^1^ = I c-W'I'l J (21) 

This corresponds to the propagation of the excitation with negative energy (hole). For the elements of the 
Green function we have at z < 0: 

Gu(iz,x) = 2ir 2 i J pdp,J (\p\x)cxp(-v\z\\p\ J ) 

Gi 2 (iz,x) = -2n 2 i J pdpi J Jj{\p\x)exp(-v\z\\p\ J ) (22) 
Next, we rotate back the Green function to real time t — iz, and get: 

Gn(t,x) = -2sign(i) 7r 2 i / pdp J (\p\x) exp(-iv\t\\p\ J ) 



Gi 2 (t,x) = -2tt 2 i J pdpi J Jj(\p\x)exp(-iv\t\\p\ J ) (23) 
Below we shall describe several particular cases. 



3.3. The case J =1 

In this subsection we consider the (emergent) relativistic invariant interaction term corresponding to 
O a — 03 that gives 



F (u;) ~ 2Re / rfte— « (24) 



1 

{v 2 t 2 ~x 2 {t)) 

and the case Oa — Va&b corresponding to Raman scattering. In the latter case 



W ~2Rey o (25) 

Substituting here the trajectory Eq. ([5|) with x(0) — we get the singularity. Therefore, we consider the 
motion along the following trajectory, which experiences the thermal- like radiation |49| . This is 

x[t] = tuoth— (26) 

Uq 

with the limiting velocity which is smaller than the effective speed of light, uq < v. The correlation function 
for the electron - hole pairs has the poles at t — and at 

t = ±^[i7r/2±(l/2)log^^] (27) 

a V — Uq 

(compare this with Eqs. A. 2, A. 7 in [49].) Therefore, both expressions for the click rate Eq. (|2~4"|) and Eq. 
(|25| contain Boltzmann exponent. Expressions Eq. (f24|) . (|25|) are divergent at t —> 0. As it was explained 
at the end of section 13.11 we apply the conventional regularization of the Wightman functions shifting the 
argument t to t — ie with small e (see also [HI). In addition, the real part of 2 is equal to . The 
contour is closed either in the lower half of the complex plane or in the upper half of the complex plane 
depending on the sign of u>. 

For the case of the emission of the electron - hole pair the dominant nonexponential term ~ w 3 is given 
by the pole at t = 0. For the case of the absorbtion of the pair the pole at t = does not enter the expression 
for the click rate, and the term corresponding to the second pole dominates. The corresponding tunneling 
exponent is given by 

Fo(w)~exp(~), lu>0 (28) 

with the temperature T — 2 ^ Uo • This is similar to Unruh temperature but with the limiting speed uo on the 
trajectory instead of the effective "speed of light" v. 

3.4- J — 1- true relativistic case 

In this subsection we consider the situation that cannot be realized in real graphene. However, it resembles 
the true relativistic case and, therefore, deserves the consideration. Namely, we suppose, that the internal 
degrees of freedom of the detector are invariant under the Lorentz symmetry. This occurs in true relativistic 
2 + 1 model. 

We consider the hyperbolic motion 

2 2 

t(s) = - sh {as/v), x(s) = — ch (as/v) - — (29) 
a a a 

with the linear acceleration a. Instead of Eq. (jT2"j) we use 

POO 

F t (w) = 2Re / dS[f]e- iuS W W{t , t - /) (30) 
Jo 

where S is the internal time of the detector, and 

/•OO 

FoH~Re / dS[t]e-' lujS ^Ti-g(-t,-x)G(t,x) (31) 
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instead of Eq. ffHl) . We have 

fOO 



F (ui) ~ Re / ^e--——!—-—- / dse — ~ ,:<•>) 



1L1 

2c 



The regularization described in section [3~T1 applied here would lead to the values of F t {ui) depending on t. 
Fortunately, for t = the answer is obtained that coincides with the result obtained in Refs. [45j, [23 1, when 
the finite small size of the detector is taken into account. We suggest the following resolution of this puzzle. 
Since the interaction between the detector and the electron - hole pairs in this case is completely Lorentz 
invariant we should consider the Wightman functions that describe emission/absorbtion in the reference 
frame moving together with the detector. The arguments of section f3. II then lead us to the conclusion that 
the shift is to be applied to variable s (time in the moving reference frame) instead of t. Then the regularized 
Wightman function W(t[s\ — ie]. t[$ 2 + ic]) coincides with the Wightman function regularized as in Eq. (20) 

of Hi. 

Thus, we assume that the integration contour is shifted to s — ie with small e. Then the integral is equal 
to the sum over poles: 



F (w) 



V 



„ lyl i. 3 M 

e T u -1 



(33) 



Here Tjj — a/(2irv) is the Unruh temperature. Thus the ratio of the absorption probability to the 
emission probability is ~ e -^/Tu _ 

3.5. The Dirac- Galilean case J = 2, the interaction term corresponds to Raman scattering 

In this subsection we consider again the interaction term Vinteraction ~ ip + a a a b ip that corresponds to 
Raman scattering. For the parabolic trajectory Eq. ([5]) the problem can be solved exactly. Here 

Qu(iz,x) = -sign(z) 2tt 2 —\- -e~ T ^ x 

Zv\z\ 

G 12 {iz 7 x) = -2tt 2 , ? , „ (-4:v\z\ + e~~^ x2 x 2 + 4e~^i x2 v\z\) (34) 
2v\z\x z 

and 

4 

W{iz) = g n (iz,x)g n (-iz,-x)~ ^—e-£-y (35) 

v z z* 

Next, we substitute z — it, and the trajectory x[t] from Eq. © with x[0] = 0: 

w(t) ~ -4r^ vFH * (36) 

The click rate (response function) has the form 

F (w) ~ - f°° dse- ius l T u ) -I e ^ 3 . (37) 

J— oo s 

(2) 

Here T^y is the Unruh temperature for vacua with J = 2, it is determined by the generalized acceleration 
in Ea.(|TD|l and is given by Eq. Q: 

'H 2 F 2 \ i/3 



T (2) _ 
1 u — 



(in this expression we have restored the Planck constant). 

Expression Eq. ([37| is linearly divergent at s — > 0. The regularization described in section [3~T1 effectively 
results in the shift of the integration contour to s — ie with small e. This is equivalent to the following regular 
expression for the click rate: 

FoM ~ .(l-signH)W r&-^(^ 3 -i). (39) 

J. jj J — OO 
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In principle, the integral can be taken and expressed through the generalized hypergeometric functions 
[53j . However, we do not present here the corresponding lengthy result since the regular expression of Eq. 
([39| itself can be used effectively for the calculation of the values of the distribution and for the drawing of 
the plot. This distribution is highly asymmetric. For uj — > — oo it tends to infinity, while for uj — > +00 it 
tends to zero. 

Alternatively, this integral can be taken using the stationary phase approximation. At uj > and 
uj ^> T^p there are the stationary points s — ± / —^37- The asymptotic expansion at large uj may be 

I — , V 3T u 

obtained similar to that of the Airy function 48] : 

,(2) 



3T{f > \ 5/4 



2uj \3/2 



F (u>) ~ ( ^TT ) cos ( T^T ) . " > 4 2) > (40) 



The derivation of this result admits the following interpretation. If we disregard the pseudospin structure, 
in the reference frame moving together with the detector the one - particle hamiltonian has the form % = 
—2vFp x t + vp x + vpy. It depends on time t, and the one - particle problem becomes non - stationary. In the 
above expressions the motion starts at t = and ends at t = t. The momentum p does not depend on time 
but it depends on t at the end of the time interval. The "classical" motion corresponds to the extremum of 
J Q W(t)dt = —vFp x t 2 + vp x t + vp 2 t with respect to p x ,p y , and gives p c i. y = 0,p c i, x = ^t. The corresponding 

energy as a function of time satisfies J* E(t)dt — — ^-t 3 . In semiclassical approximation the wave function 

of the electrons/holes contains the factor exp(— i J* E(t)dt) ~ exp(i^-t 3 ). Because of the interaction with 
the detector (the interaction potential is given by Eq. (fTTj) ]) the hole may be transformed to electron at any 
time. Therefore, the process is allowed, when first, the pair is created, and second, it is annihilated. The 
amplitude is proportional to exp(— 2i J Q E(t)dt) ~ exp(i : ^-t 3 ) . The amplitude of the process with the two 
changes of the energy of the detector by ±w is given by the Fourier transformation 

A~ J dt(,..)exp(-iut-2i J E{i)dt) = J dt(...) exp(-iujt + i^-t 3 ) (41) 

By unitarity the total probability that the detector emits/absorbs the pair is proportional to the imaginary 
part of A. 

Signs in Eq. ([41]) are important. This is also important that positive values of u> correspond to the 
absorption of the electron - hole pairs. Those signs can be checked as follows. Eq. (14 1 p is similar to the 
quantum mechanical expression for the amplitude of the excitation of the quantum level with energy 2E 
by the absorption of the light with frequency — uj > 0. In the simplest case when E is positive and does 
not depend on time the integral over t gives the delta - function of (u + 2E) . This delta - function may be 
nonzero for the excitation of the quantum level (w < 0). It vanishes for the inverse process (u> > 0) because 
the system is assumed to be in the ground state. 

At negative values of oj corresponding to the emission of pairs and at \u>\ 3> T^y the dominant term 
corresponds to the Loran expansion at s = 0. This gives the power - like behavior: 

FoH ~ 2 ^^i, uj<0. (42) 

Equations (1371) and (|40[) are similar to equations which determine thermal effects in systems at quantum 
criticality, where the thermodynamic and kinetic parameters depend on frequency and temperature through 
their ratio ui/T. This generalizes the Hawking and Unruh effects to the vacua with anisotropic scaling. This 
consideration clearly demonstrates an intimate relation between Lorentz invariance of the field theory and 
Boltzmann distribution in statistical physics. For the theory at J — I with emergent Lorentz invariance the 
trajectory of the detector Eq. 0$ being substituted to the expression for the click would give the singularity. 
Therefore, we considered the alternative trajectory Eq. (|2"B1) that leads to the Boltzmann distribution of 
the pairs absorbed by the detector Eq. (|28|) . The true Lorentz invariant theory leads to the Boltzmann 
distribution as well. For non-Lorentz invariant theory with J = 2, there is no Boltzmann distribution for 
absorbed particles. Instead, we have the law of Eq. (|4"0j). 
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3.6. Arbitrary J > 2, the semiclassical consideration 

For arbitrary J > 2 we have (in the case of Raman scattering with Oa — o~ a <7b)'. 

Fo(w)~(7o(w) ~ / dte- iut g 11 (t,x)G 11 (-t,-x) 

p oo — lO poc poo 

= 4?r 4 / dt pdp p'dp' J {px[t}) J (p'x[t])ex.p{-iojt - ivt(p J + {p'] J ), 

JO-iO -JO JO 

x[t] = vF J ~H J (43) 

(For the definition of cr(ui) see subsection 13.71 ) 

Already for J = 3 the elements of the two - point Green function are represented as a rather compli- 
cated composition of generalized hypergeometric functions. Therefore, for J > 3 we apply the semiclassical 
approximation that works at |w| S> Tu. One of the advantages of the semiclassical consideration is that the 
pseudospin degrees of freedom of the electrons and holes are neglected. Therefore, the practical choice of 
operators Oa entering the interaction term Eq. pip does not matter. 

Let us consider the trajectory Eq. §5§ for J > 2. The dimensionality analysis demonstrates that the 
general form of the click rate is 

/oo ... 
dse-^ s ' T v gj {s) (44) 
-oo 

with some function gj(t), and given by Eq.([7]): 

T^ J) = (F J vyK J+ V = a/v F . (45) 
At w > Tjj (absorption) in semiclassical approximation we have 

G(t,x[t}) - E ^exp(*i^,FV+V-W), (46) 

N=0,l;K=0..J-2 ' 

where Cnk are constants, and 

/•oo 

F (u>)~Re / dte- iu,t g 2 {t : x[t]) (47) 



There are several stationary points in the complex plane. The dominant contribution corresponds to the 
terms with N = K = in Eq. (|46l) . Stationary phase approximation gives the following asymptotic behavior 
(the consideration is similar to that of the Airy function (48j): 



const 



1 u 



jJ/(J-i) 

P = [2(J 2 -l)]i/./(J+l) (48) 

This result corresponds to the absorption of the electron - hole pairs by the detector. Emission of the 
pairs by the detector corresponds to the power - like behavior that can be extracted from the Loran expansion 
of the integrand in Eq. (|47|) . Since the integration contour should be shifted to t — iO we can calculate the 
residue of the integrand at t = 0. The result gives the power - like dependence on uj: 



F„(o;) ~ (Jgy) , M»4 J) > w<0. (49) 

Applying the same technique to Eq. (|32p we can obtain the dominant power - like term for the emission 
of the pairs in the relativistic case ~ |u;| 3 that can also be seen in the exact result Eq. (|33|) . This behavior 
is in some sense similar to that of considered in (Hoj and in Sec. 31.4 of Q, where the power-law asymptote 
for radiation from rotating black hole is obtained by tunneling method. 
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3. 7. Discussion of possible experiments with the accelerated detector of electron - hole pairs 

Experiments for the case of "accelerating detector" look difficult but not impossible: one should create 
a laser spot moving along the graphene surface with acceleration and measure the Raman scattering from 
such a spot. The emission (adsorption) of electron-hole pairs by the moving spot corresponds to Stokes 
(anti-Stokes) Raman scattering, with the lower (higher) frequency of the scattered light in comparison with 
the incident one. In order to check the distributions obtained in this section this spot should move along the 
corresponding trajectories given by Eq. ([5]). One should expect a huge asymmetry, with the oscillating factor 
for the dependence on to for the absorption versus the power-law for the emission. Unruh effect corresponds 
to the absorption of the electron - hole pairs by the accelerated detector, or, better to say, the Unruh effect 
probability is the ratio of the absorption probability to the emission probability. In order to model the Unruh 
detector the absorbtion of the incoming photon and emission of the secondary photons should occur in the 
small vicinity of the spot and without the essential delay. Otherwise the spot cannot be considered as the 
point - like detector moving along the given trajectory. 

The detailed consideration of this process, and the questions related to the delays and the correlation 
length of Raman scattering from the spot is out of the scope of the present paper. Nevertheless, in order 
to understand, what plays the role of the operator Oa of Eq. (fTTj) in this case, let us consider the simplest 
diagram that describes the Raman scattering. The incoming photon is transformed to the electron - hole 
pair. Next, the electron and the hole scatter elastically on the atoms of the crystal lattice. During this 
scattering process the momentum of the incoming photon (minus the momentum of the outgoing photon) 
is absorbed by the crystal lattice. After that the scattered electron emits the outgoing photorQ The initial 
state is vacuum, the final state includes the created electron - hole pair. 

For the cross section of the Raman scattering from the moving spot we have 



Here the sum is over the polarizations a, b = 1,2 of incoming and outgoing photons. The so - called click 
rate is given by the derivative of a with respect to T in the integral over t: 



We come to Eq. (|T2"j) calculated for the interaction term Eq. ffTT]) with Oa — o- a o~b- Eq. ([5"Tj) was derived 
for the emission of electron - hole pairs. This corresponds to negative values of a;. For the spot at rest o{ui) 
vanishes for uj > since there are no free electron - hole pairs in vacuum. However, for the accelerated spot 
ct(lu) already does not vanish for positive lo. This means that the energy of the scattered photon may be 
larger than the energy of the incoming photon. This can be interpreted as the detection of the electron - 
hole pairs by the Unruh detector. 

It is worthwhile to note that the moving laser beams were successfully used to study quantum phenomena 
in ultracold gases [Hl|. The laser spot moving along the circular trajectory in BEC has been suggested (H^j 
to study the analog of Zel'dovich-Starobinsky effect which is similar to the effects discussed here. 

4. The system in the presence of electric field 

In the presence of external electric field electron - hole pairs are created due to the Schwinger mechanism. 
These pairs may be detected in some way by the detector at rest. This situation is somehow similar to that 
of the accelerated detector considered in the previous section. Namely, free electrons and holes would move 
with the acceleration (in opposite directions) . In some sense we may speak on the "accelerated vacuum" in 
this situation. 



2 Actually, the hole may emit the photon, and this may happen before the scattering on the crystal lattice. Also, the photon 
may first scatter elastically on the lattice, and then be transformed to the electron - hole pair. But the description of the 
corresponding processes will lead to the same final answer for the operator Oa- 




dt / dt'e-^-^ ^(0|V + (x[t']) ( T^ a V(a:M)V' + (xM)Ta^^(^M)|0) 



(50) 




(51) 
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Distribution of electron-hole pairs as a function of energy 

The one-particle Hamiltonian for 2D massive fermions with anisotropic scaling E 2 = v 2 (p 2 + p 2 ) J in 
electric field £ has the form 

H= j J \ (52) 

V[jJ X + iPyj £X J 

Here pj = —idj] we use the units U = e = 1. 

Integration over the classically forbidden region gives us the pair production probability (probability of 
the Schwinger effect) [§] 



hoi 2 



exp {^ a jP'l +1 \. 



The energy of the electron-hole pair is u> — 2v(p 2 +Py) J ^ 2 - The distribution of the created pairs as a 
function of energy has the form: 

/H = J dp x dp y e- a TPi +1 6(u;-2v(p 2 x +pl) J / 2 ) 

(#) ^ex P (-a((-) -a 2 ) ) 



( V H 77 \ 2 / J \ 



1/2 



dx- 



(J+l)/.// x(J+l)/2 

exp(-a(^) (l~x 2 



1/2 



dz exp(-a(| ; ) (J+1)/J z)z 1 /(^i)-i( 1 _ z2 /(./ + i)) 1/2 - 1 (54 ) 



Let us introduce the Unruh temperature T = in correspondence with Eq. ([7|) : 

T« = (f ^) 1/(J+1) (55) 

(we use here £ instead of F.) 



4--2. Particular cases 
For J — 1 we arrive at 



/(«) ~ J'dzexp "M '^-'(l-z)" 2 1 



This distribution was also derived in Ref. 54j . where the spontaneous recombination of electrons and 
holes created due to the Schwinger effect in monolayer graphene was considered. 

For J > 2 the function f(uj) is expressed through the generalized hypergeometric functions (for their 

definition see |53|). In general case w 2 /(w) is the analytical function of ( = a( ^ J : 



^ 2 /M~C r^exp^C^ 17 ^^-^!-^/^ 1 )) 172 1 (57) 
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This integral is convergent for any values of J. It can be represented as a composition of the generalized 
hypergeometric functions. For example, for J = 2: 



w 2 /( w ) 



4 1*. 

3' ' 3' 

5 7 3 3 

6 ' 6 ' 2 ' 2 



+ |C2^3 



1 5 
,6' 6' 

- - 1 

3' 3' x 



In this expression the functions p-Fg at small values of the argument are given by 53 1 



(X 1 , . . . , dp , £ 

h,...,b q 



n=0 



T^n 7T1 T' W»=¥ + l'> + 11 - 1 

(oi) n ...(& g )„ n! 



(58) 



(59) 



4.3. Asymptotic expansions 

The expansion of Eq. (|59p works in a vicinity of z — and does not work for z — > oo. In order to evaluate 
the distributions at £ — > oo we use integral representation Eq. ([57} . We obtain: 



/(") 



C^ra/cj + i))^!^ 



w > r. 



!7 



(60) 



One can see that the dominant term is power - like. Among the other terms there are also the subdominant 
terms that contain the tunneling exponent (this may be obtained as an expansion near z = 1 in the 
integral of Eq. (|57p). It is instructive to compare this with the behavior of the distributions of electron - hole 
pairs detected/emitted by the moving Unruh detector. Eqs. (|57l) . (|60)) are to be compared with Eq (|48|) . In 

( u \ 1/J+1 

both cases there is the dependence on £ = ( \ . This indicates that the two phenomena may indeed 

be similar, but there is no direct correspondence. In the distribution of the electron - hole pairs created in 
external Electric field the leading term is power - like. This is somewhat similar to the distribution of pairs 
absorbed by the accelerated detector Eq. P5|) . However, the powers are different, and in the latter case the 
distribution contains also the oscillating factor. 

4-4- /( w ) expressed through the Green functions 
Expression Eq. (|54j) can be rewritten as 



/(") 



aire' 



d Px d Py e- a ^ p 'y +1 -' T2v ^ +p ^ J/2 ~ / dre^ T W(T) 



Recall that the expression Q = ~afp^ +1 has its origin at 



Q = 2 Im / Vdx, 



where 



V 



\J([E-£x)/v) 2/J -pi 



(61) 



(62) 



(63) 



is the classical momentum. This, in turn, may be obtained as a semi-classical approximation to the repre- 
sentation 



(64) 



f(w) = Re / dTe™ T (ip(0,0)1)(0,0)ip(T,x)ip(T,x))dx 
- Re J dTe tulTr Trg £ (T,x)g £ {-T,-x)dx 



Here fermion Green functions in external electric field are present. This expression is to be compared with 
Eq. (H5). 
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4-.5. Discussion of possible experiments with the "accelerated vacuum" of graphene 

Relevant energies of electron excitations for this case lie in the infrared, which allows us to suggest for the 
detection an efficient optical tool, namely, the electron Raman scattering (see, e.g., 57|). One should probe 
the dependence of the Raman spectra on electric field. This will correspond to our case of "accelerating 
vacuum" . Another way to probe the distribution of the electron - hole pairs on energy is to consider the 
spontaneous recombination of electrons and holes forming the electron - hole plasma [54| which may be in 
principle probed via luminescence. 

Of course, estimation of the Unruh temperature is crucially important for the choice of the proper exper- 
imental technique. First, notice that graphene is very robust material with respect to electric breakdown; 
the current densities such as (3 — 4) x 10 4 A/m are reachable (for recent references, see [EH (56|) which 
corresponds to the electric fields £ c m 10 s — 10 9 V/m. According to Eq. (1551) . for the case of bilayer graphene 
it corresponds to Unruh temperature as large as 0.5 eV ~ 6000 K. So strong electric fields produce a strong 
heating [56J and this, of course, can be a serious problem but even for the orders of magnitude weaker fields 
the Unruh temperature can be still quite high (it is about 10 meV for £ sa 3 x 10 6 V/m). Relevant energies 
of electron excitations for this case lie in the infrared, which allows us to suggest for the detection an efficient 
optical tool, namely, the electron Raman scattering (see, e.g., [HtJ)- 

One should keep in mind, however, that below the energy of the order of 10 meV the effects of trigonal 
warping and electron-electron interactions become important 0, [3^ | with the reconstruction of parabolic 
touching point to several conical touching points at low energies. This means that for the electric fields 
£ ~<C 10 6 V/m) we will probe these conical points rather than the parabolic point, and the Unruh spectra 
should be similar to those for the single-layer graphene. However, the Fermi velocity for these new conical 
points is much smaller than in the single-layer graphene (1.4 x 10 5 m/s and 9 x 10 5 m/s, respectively 
[391 ]) which will lead to an essential increase of the Unruh temperature (inversely proportional to the Fermi 
velocity), the effect which looks very interesting by itself. 

Alternative types of detection of the created electron-hole pairs and their spectroscopy are emission 
Mossbauer spectroscopy (one should put on graphene, e.g., radioactive 57 Co) [HI] or electron spin resonance 
[59] ] of magnetic adatoms on graphene. In both cases the spectra are modulated by the fluctuating electron 
density. However, optical methods (Raman spectroscopy) look the most promising and realistic. 



5. Massive ID fermions 

In this section we consider the dynamics of the one - dimensional system with the Hamiltonian given by 
Eq. (fTJ), where p y is fixed and is equal to M playing the role of effective mass of the fcrmionic excitations. 
The multilayer graphene itself (at ABC stacking) can be considered as a collection of such ID systems 
corresponding to different values of M . The one - dimensional system considered in this section is naturally 
realized in graphene nanoribbons that is the thin strip of the multilayer graphene. The values p v are quantized 
there [9j, and the subsystems with different values of M — p y can be distinguished from each other (see, e.g., 
the discussion of conductance quantization in Section 5.5 of the book Q). 

5.1. Miiller temperature for 1+1 massive fermions in electric field 

In the presence of constant external electric field we deal with the "accelerated vacuum" . We consider 
the Hamiltonian Eq. (JT]) in the case when p y is fixed and equal to M. Then p — p x + iM, and one-particle 
Hamiltonian for ID massive fermions with anisotropic scaling E 2 — v 2 (p 2 + M 2 )' 1 / 2 in electric field £ has 
the form 

»=l , v { ' I 

v[px + iM £x 



Integration over the classically forbidden region gives us the pair production probability [8| 

l9 -2JB\ § ,i)vM J + 1 /£ , 

\vo\ 2 = e V 2 ) (66) 
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With the definition of generalized acceleration a = v F £/M of Eq.© this expression resembles the thermal 
distribution. Namely, we may represent 



vM J 
Tm 

_ lt , p £ 2J n ^3 J 



|?7o| 2 = exp 

T u =l- L ha/v F , a = v F jj=v F ^ , 7 =riB(J,|j. (67) 

with the analogue of Miiller temperature [2(| Tm- (In the last expression we restore Plank constant h.) 
In the relativistic systems (J = 1) the Miiller temperature Tm = ha/irc is twice the Unruh temperature 
Ttj = Ua/2'nc 17j, see e.g. recent paper [2l|. The only difference from the relativistic system is that in 
anisotropic scaling system the pre-factor depends on J. 

5.2. Unruh effect at \ui\ Tjj 

In this subsection we consider the situation, when the detector moves in the given ID system along the 
trajectory given by Eq. (j4}. With p = p x + iM one has the Green's function of ID quasiparticles 



G{it,0) = 2m I dp x e 



JPxX-v\t\\p\ J 



1 1 si g n wfe) J 

U gn(<) te) J 1 



Gii (it) G 12 (it) 
G 2i (it) G u (it) 



Gii (it) = (2tt)iJ d Px exp(i Px x-v\t\(pl + M 2 ) J ' 2 ) 

G i2 (it) = sign(t) (27r)i j dp x expfex - v\t\(p 2 x + M 2 ) J ' 2 ) ( ^ X ~^p )1/2 ^ ' 
G 2 i(it) = sign(t) (27r)i J dp x expfex - v\t\{p 2 x + M 2 ) J / 2 ) ( j^fj^IJi 
We evaluate the Green function in the semi - classical approximation. Then 

G(t,x[t]) ~ J dp x e ip * x - iv WW J (69) 

Stationary phase conditions give 

x = v\t\Jp x (p 2 x + M 2 ) J / 2 - 1 (70) 

Below we consider the two particular cases (J = 1,2) of Unruh radiation of massive particles with the 
energy near to the threshold 2E(p = 0) = 2vM J . 

1. J= 1. 

In this case we have 



e .>r/4 

As a result 



CHI. .r[/-]) - ——^(-ivM^jl- (^) 2 \t\) (71) 



F (lj)~RcJ dt-^exp(~iujt - 2ivMy 1 - (^)^) ~ ^ dt-^exp^-iwt - 2ivM^ 1 - (^) 2 *) 

(72) 

Here the trajectory of the detector x(f) is given by Eq. (Q}: 

x(t)=v(t 2 + ^) 1/2 -v- (73) 
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We also use the following parametrization of this trajectory: 



We come to 



t(s) = - sh {as/v), x(s) = — ch (as/v) - — (74) 
a a a 

f°° 1 / iv r as i \ 

F (w)~i / ds— — rrexp u sh (as/v) + 4vMsh — (75) 

J_ 00 th (as/v) V a L 2viJ 



At w ~ 2wM 3> T[/ = -f^j we can apply the stationary phase approximation. For < cj < 2vM and 

2«]Er << 1 there are no stationary points at real s but there is the stationary point at imaginary 
values of s in the lower half of the complex plane. Therefore, the absorbtion of the electron - hole pairs 
is described by the tunneling exponent. One can define the effective chemical potential /i e // — 2i>M 
and the effective Unruh temperature T\j — = ^p-- Then 

F (w)~e V r ^ / , \u - 2vM\ < 2<uM, w > o/u, w < 2uM (76) 

It is worth mentioning that in this case the Unruh temperature is not given by Eq. ([7]). Instead we 
have the relation that involves both mass parameter M and the acceleration a. For u> > 2vM there 
is the stationary point for the real value of s that results in the appearance of the oscillating factor in 
F (uj): 

1/2 



F (w) ~ sin (AVE r Me// ) ) , |w - 2t>M| < 2«M, w > a/v, cj > 2vM (77) 



J = 2. 

In this case 



As a result 



0(t,a:[f]) ~ ^exp[-v\t\(Ft-M)(Ft + M)) (78) 

f (w) ~ / di— exp(-iu;i + -vt(Ft - M)(Ft + M)) (79) 
./-^ V 2 / 



Next, we use the stationary phase approximation. It is valid at |w| ^> Tjj = (vF 2 ) 1 ^ 3 . Positive values 
of uj correspond to the absorbtion of the electron - hole pairs. The Unruh temperature for this case 
Tjj ~ (vF 2 ) 1 / 3 is equal to that of defined by Eq. (|7J| with the generalized acceleration of Eq. (fTU|) as 
expected for the J — 2 case, when the generalized acceleration does not depend on M. We denote 
fj, e ff — 2vM 2 , and get similar to the case of the Airy function [481 ] : 

• . , /2v / 6/w + ixeff \3/ 2 \ \lo-2vM 2 \ , . 

5.5. Suggested experiments 

To probe the distributions obtained in this section we suggest the same methods as suggested in Sections 
III and IV but for the graphene nanoribbons. Namely, Raman scattering and spontaneous recombination 
can be used to probe the spectrum of electrons and holes created during the Schwinger process. The Raman 
scattering from the moving spot can be used to probe the Unruh effect itself. 



6. Conclusions 

Vacua with anisotropic scaling, such as r — > br , t — > b z t, represent new topological classes of quantum 
vacua, where the scale parameter z is determined by the topological charge of the quantum vacuum. The 
relativistic invariance emerges only in one of the the universality classes of vacua, where topology supports 
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the isotropic scaling z — 1. Anisotropic scaling naturally emerges in many condensed matter systems, 
including multilayer graphene. That is why it is instructive to consider general properties of such quantum 
vacua. Here we discussed the quantum effects, such as Unruh radiation and Schwinger pair production, in 
vacua which belong to universality classes with topologically protected anisotropic scaling. Though these 
effects do not possess thermal properties relevant for the relativistic class of quantum vacua with z = 1, 
some features are proved to be universal. Both effects, Unruh radiation and Schwinger pair production, are 
characterized by the modification of correspondingly Unruh and Miiller temperatures. This characteristic 
temperature is proportional to the acceleration properly generalized to systems with the anisotropic scaling 
in Eqs. ^ and (|10[) . As a rule such temperature does not describe any thermal distribution of radiated 
excitations. Instead it acts as the temperature, which enters the thermodynamic and kinetic parameters of 
the many-body system in the vicinity of the quantum phase transition, where the vacuum experiences the 
properties of quantum criticality. 
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